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We report on the experimental observation of beam-like neutron emission with peak flux of the 
order of 10® n/sr, from light nuclei reactions in a pitcher-catcher scenario, by employing MeV ions 
driven by high power laser. The spatial profile of the neutron beam, fully captured for the first 
time by employing a CR39 nuclear track detector, shows a FWHM divergence angle of ~ 70°, with 
a peak flux nearly an order of magnitude higher than the isotropic component elsewhere. The 
observed beamed flux of neutrons is highly favourable for a wide range of applications, and indeed 
for further transport and moderation to thermal energies. A systematic study employing various 
combinations of pitcher-catcher materials indicates the dominant reactions being d(p, n+p)^H and 
d(d,n)^He. Albeit insufficient cross-section data are available for modelling, the observed anisotropy 
in the neutrons’ spatial and spectral profiles are most likely related to the directionality and high 
energy of the projectile ions. 


Neutrons provide many opportunities for probing of 
materials in ways which charged particles and ionizing ra¬ 
diation cannot. In this context an ultrashort, directional 
burst of neutrons with high flux in the MeV range would 
have wide ranging applications. Exciting opportunities 
for ultrafast studies lie in the area of materials for fusion 
energy research due to growing interest in understanding 
neutron-induced damage at the atomic scale [l[ . Further¬ 
more, a compact source of pulsed, MeV neutrons would 
provide novel capabilities for interrogation of large camo 
containers by fast neutron radiography techniques [I, [H , 
where the nature and location of the threat can be identi¬ 
fied by simultaneously measuring scattered neutrons and 
time of flight of the induced gamma radiation. 

Significant attention has been paid recently to laser 
driven sources capable of producing short neutron bursts, 
and having potential advantages in terms of cost reduc¬ 
tion and compactness, reduction of radioactive pollu¬ 
tion and ability of radiation confinement by close-coupled 
experiments. Although a different approach (high en¬ 
ergy deuteron-breakup) has recently been reported Q, 
the most established route to create a laser based neu¬ 
tron source is by employing laser accelerated ions in ei¬ 
ther fusion or spallation reactions. Since spallation of 
heavy atoms requires high energy projectile ions, reac¬ 
tions based on low atomic mass nuclei, such as protons, 
deuterons, lithium etc., are particularly relevant. The 
neutron yield from nuclear reactions scales with the prod¬ 
uct of the densities of the interacting species and the 
cross-section a, which for most common reactions reaches 
high values at ^MeV centre-of-mass energy. Producing 



FIG. 1. The experimental set-up showing the in-chamber 
diagnostics and arrangements of nToF detectors around the 
chamber. 


high fluxes of MeV ions using intense lasers is currently a 
very active area of research. Where a number of emerging 
ion acceleration mechanisms, such as radiation pressure 
acceleration (RPA) @ and breakout afterburner (BOA) 
0, hold the promise for producing higher energy ions 
with higher efficiency, target normal sheath acceleration 
(TNSA) 0 is a well-established and robust mechanism 
which produces MeV ions with high flux and narrow 
divergence. Such beams can be readily deployed in a 
pitcher-catcher setting for neutron generation via nuclear 
reactions. 

In addition to the benefit of high reaction cross- 
sections, anisotropy in neutron emission is another facet 
of beam-nuclear reactions. Simulations [H show that us¬ 
ing several MeV ions in the above mentioned reaction 
involving low Z materials yields a neutron flux strongly 
peaked along the ion beam forward direction, and that 
the anisotropy grows further with increase in ion beam 
energy. A beamed neutron source is highly favourable 
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FIG. 2. The neutron beam can be clearly seen in the lower right of the CR39 slide in (a). A close-up of the pits created on the 
CR39 can be seen in (b), obtined by using a transmission microscope with lOx zoom, which can be compared with the (d), in 
which there is no catcher target to convert the ions into neutrons. The laser axis was 1 cm up from the bottom of the CR39, 
and 1 cm in from the right hand edge, (e) shows the horizontal and vertical beam profiles from the scanned CR39 slide, along 
with the neutron flux (between 2.5-10 MeV) obtained from the nToF and bubble detector data from the same shot. 


not only for the aforementioned range of applications, 
but also for its further transport and an efficient moder¬ 
ation to thermal energies required for another range of 
applications including for instance, Boron Neutron 
Capture Therapy Hi- Whilst some degree of neutron 
beam anisotropy has been recently reported at low neu¬ 
tron fluxes [ll|,ll2|, here we show, for the first time, direct 
imaging of a true beam-like neutron emission with peak 
flux of the order of 10^ n/sr. In particular, a neutron 
beam with a FWHM of (70 ± 10)° and peak flux of (5 
± 2) X 10^ n/sr was captured in nuclear track detector 
kept in the proximity to the source. The neutron beam 
was produced by d(p, n+p)^H and d(d,n)^He reactions 
driven in a beam-catcher scenario, by employing protons 
and deuterons from thin deuterated plastic (CD) foils ir¬ 
radiated by a sub-petawatt laser. 

The experiment was performed using the petawatt arm 
of the Vulcan laser at the Central Laser Facility of STFC, 
UK. The linearly polarised laser pulse of 1053 nm wave¬ 
length, delivering ^200 J on target after reflection from a 
plasma mirror, was focused down to a focal spot of ^6 jam 
FWHM, providing peak intensity on the target 3x10^0 
Wcm'^. Various targets were irradiated by the laser in 
order to generate energetic ions via the TNSA mecha¬ 
nism, namely gold foils (10 fim thick), 98% deuterated 
polyethylene ((72^4)^ (henceforth called CD) foils with 
and without few microns thick A1 foil at the rear side. 
The ion beams were diagnosed in the earlier part of the 
experiment using a In^h resolution Thomson parabola 
spectrometer (TPS) [13 with image plate detectors look¬ 
ing at the target normal direction, as shown schemati¬ 
cally in Fig. [TJ Since traces for ions with the same charge 
to mass ratio overlap in the TPS, we implemented the dif¬ 


ferential filtering technique described by Alejo et al Q 
in order to extract the deuteron spectra from the diag¬ 
nostic. ^2 mm thick solid blocks of CD and graphite 
were placed ^ 3 mm behind the pitcher target (hence¬ 
forth called catcher) in order to generate neutrons from 
nuclear reaction. The transverse size of the catcher was 
large enough to capture the entire ion beam. A full suite 
of neutron diagnostics was deployed in order to diagnose 
the spatial and spectral profiles of the neutrons generated 
in different shots. In order to capture the ffux profile of 
the emitted neutrons over a large solid angle, and due 
to its low detection efficiency, of the order of 10“"^, the 
CR39 nuclear track detector was placed in close proxim¬ 
ity (7 mm) to the catcher. The CR39 was shielded by 4.5 
mm thick lead in order to stop the high energy protons 
(upto 50 MeV M) produced at the pitcher target from 
reaching the CR39. The absolute neutron flux was ob¬ 
tained from the CR39 by using the etching method and 
calibration given by Frenje et Absolutely cali¬ 

brated Bubble Detector Spectrometers (BDS) [l7| and 
nuclear activation diagnostics [HI were used, behind the 
CR39 along the ion beam forward direction, at a dis¬ 
tance of 10 cm and 50 cm respectively from the pitcher 
target. Since the BDS provides absolute neutron flux 
in six discrete energy intervals within the energy range 
from 0.01 MeV to 20 MeV, it is possible to ascertain the 
flux of MeV neutrons generated in the catcher by dis¬ 
counting the large signal produced by the lower energy, 
scattered neutrons bouncing within the target area and 
hitting the detector several times over a period of time. 
Using 1 mm thick pure indium foils and measuring the 
decay of ^^^In and ^^^In, the neutron flux was estimated 
within two energy intervals, namely 0-3 MeV and 0.7-15 
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FIG. 3. (a) and (b) show the neutron flux (from 2.5 MeV to 
10 MeV) along (a) the laser axis (average of the nToF detec¬ 
tors at 2° and 8°) and (b) off laser axis (average of the nToF 
detectors at 77°, 139° and 148°) obtained for different pitcher 
and catcher materials and diagnostics. Where the color of the 
data points represent the pitcher and catcher material com¬ 
binations, as mentioned in the figure legend, ’x’, square and 
circle represent the data obtained by scintillators, activation 
and BDS respectively. 

MeV. Six neutron time-of-flight (nToF) detectors con¬ 
sisting of EJ410 plastic scintillators, optically-coupled to 
XP3330 photomultiplier tubes (PMT), were used to pro¬ 
vide spectral information at different emission angles by 
the time-of-fiight method. The detectors were shielded 
appropriately (typically by 5-10 cm of lead on every side 
and 5 cm of plastic all around except the front of the 
detector) in order to suppress saturation of signal from 
Bremsstrahlung radiation and to reduce the noise in the 
signal due to scattered low energy neutrons reaching the 
detector. The angle of observation and distance of each 
nToF detector is given in Fig. [TJ The nToF detectors 
were cross-calibrated 0 against the spectra obtained 
from the BDS, which were absolutely calibrated by the 
company [l^. 

As it can be clearly seen in Fig.[2{a), the CR39 shows a 
beamed neutron emission from the catcher along the ion 
beam forward direction. The (0,0) co-ordinate in this im¬ 
age represents the ion beam axis, which is near the bot¬ 
tom right corner of the CR39. For comparison one can 
see Fig. He) showing the CR39 from the reference shot 
taken without the catcher, while keeping everything else 
the same. Both shots had identical pitcher targets of 10 
/im CD and measured laser energy to within 1% of each 
other. As expected, the amount of neutron generated in 
the later case (due to the interaction of ions, electrons 
and gamma rays produced from the pitcher target with 
the surrounding objects, including the lead shielding of 
the CR39) was not significant and lower than the detec¬ 
tion threshold of the CR39. Neutrons are diagnosed in 
the CR39 due to the latent tracks created by the knock- 
on protons, which are revealed after etching in an alkali 
solution. The etched pits produced in the CR39 can be 
seen in the zoomed view of the CR39s, shown adjacent 
to the respective images. Using the calibration for detec¬ 
tion efficiency of CR39 given by Frenje et. al. [HI, the 
pit density across the CR39 was converted into neutron 
flux, as plotted in Fig. [2]^e). The horizontal and vertical 


FIG. 4. (a) Graph showing neutron yield from pitcher-only 
configuration (red circles) using GD targets reported in lit¬ 
erature, as cited beside the data points. The red star rep¬ 
resents the average neutron yield from 10 /xm GD target ob¬ 
tained in our experiment, as shown in Fig. (3] (b) Graph 

showing the ratio between on and off axis neutron fluxes re¬ 
ported in literature, as cited beside the data points, with 
respect to on axis flux obtained in experiments. The data 
point obtained from our experiment is labelled as [i], [ii] and 
[hi] which corresponds to different pitcher and catcher com¬ 
binations, such as GD^GD, Au/GH^GD and GD+Al^GD 
respectively. The Au/GH^GD case produced the highest de¬ 
gree of beam anisotropy, which is expected due to the low level 
of isotropic neutron flux produced from the non-deuterated 
pitcher targets, as can be seen from Fig. O 


lineouts of neutron flux across the CR39, passing through 
the co-ordinate (0,0), show an axisymmetric neutron flux 
profile with a FWHM divergence of (70 ± 10)°. The peak 
neutron flux along the ion beam axis was estimated as 
(5 ± 2) X 10^ n/sr. 

In order to identify the dominant nuclear reactions pro¬ 
ducing the beamed flux of forwardly directed neutrons, 
a systematic study was carried out by varying different 
pitcher and catcher materials. As can be seen in Fig. [3l 
the pitcher-only configuration using lO/xm CD target gen¬ 
erates a fairly isotropic neutron emission in Air with an 
average flux ^ 5 x 10^ n/sr. This corresponds to a to¬ 
tal neutron yield in excess of 10^ neutrons (> 10^ neu- 
trons/Joule), which is in line with the trend of neutron 
yield with respect to the incident laser intensity reported 
in literature (see Fig.|4l[a)). In this case the neutrons are 
generated either by the thermonuclear reactions in the 
hot dense plasma produced by the laser interaction, or 
by the fusion reaction in the target bulk driven by the 
ions accelerated at the laser front surface through the 
hole-boring mechanism [28|. 

The isotropic neutron emission from pitcher-only tar¬ 
gets, which is commonly observed in experiments as 
shown in Fig. S^b), is in stark contrast to the forwardly- 
directed, beamed neutron flux obtained in pitcher- 
catcher configuration as shown in Fig.[3]and Fig.[5][a). As 
can be seen in Fig. [4|[b), our data represents the highest 
degree of beam anisotropy observed experimentally with 
peak flux of the order of 10^ n/sr. The order of magnitude 
increase in neutron flux along the ion beam axis (which 
is same as the incident laser axis) observed in this case 
originates from several different nuclear reactions inside 
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FIG. 5. (a)-(d) show the comparison between angularly re¬ 
solved neutron spectra obtained by the six scintillator detec¬ 
tors (shown in Fig. [T]) obtained in different pitcher-catcher 
combinations as labelled on the top of the graph. The graph 
also show the data obtained from BBS and activation diagnos¬ 
tics, where available. Where BBS and activation diagnostics 
were installed along the beam axis in case of (d), as shown in 
fig-IU the BBS in (c) was installed at 145 deg from laser axis. 


the catcher. Possible reactions include the interaction 
of protons, denteron and carbon ions produced from the 
pitcher target by the TNSA mechanism with the carbon 
and deuterium atoms in the catcher. 

Different pitcher-catcher combinations were used in or¬ 
der to unfold the contributions from different potential 
reactions towards the observed beamed flux of neutrons. 
Firstly, employing graphite catcher with CD pitcher tar¬ 
get (CD^C) did not produce (in spite of shot-to-shot 
fluctuations) a significant increase in the neutron flux or 
beam anisotropy compared to pitcher-only target. This 
suggests insignificant contribution being made by the 
p+C, d+C and C+C reactions in the context of the ob¬ 
served neutron flux in the CD^CD case. Secondly, in 
order to discriminate between proton and denteron in¬ 
duced reactions in the CD targets, several shots were 
taken by using a CD catcher in front of a deuterium 
free ion source (No_d^CD) by using either lOjam gold 
or CH foils as pitcher targets (Au/CH^CD), or 10/im 
CD foil backed by few microns thick aluminium foil as 
pitcher targets (CD+Al^CD). Similar neutron flux as in 
the CD-^CD case was observed in this case (No_d->CD). 
Since the denteron and carbon ion spectra produced from 
the pitcher targets are similar in terms of number of par¬ 
ticles and beam temperature, as shown in Fig. [6fa), d+C 
reactions in case of No.d^CD can also be assumed in¬ 
significant. Therefore, the most promising reaction in this 
case would be the proton driven neutron generation via 
breakup of deuterons (d(p, n+p)^H) in the catcher. This 
is expected due to (1) large number of high energy pro¬ 
tons (^20 MeV) being produced from the pitcher tar¬ 
get by the TNSA mechanism, as shown in Fig. Ha), 
and (2) high cross-section of d(p, n-fp)^H reaction (al¬ 
though there is a limited amount of data available [1^ ) as 
shown in Fig.[6l^b). Since cleaning techniques for removal 
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FIG. 6. (a) Gomparison between on axis proton (red), 

denteron (blue) and carbon (black) ion spectra obtained from 
10 fim thick GD target, (b) Gomparison between reaction 
cross-sections for d(p, n+p)^H (red) and d(d,n)^He (black) for 
different projectile ion (protons and deuterons respectively)) 
energies, as obtained from EXFOR [2^. 


of target contaminants (such as laser ablation, resistive 
heating etc.) could not be implemented in the experi¬ 
ment due to setup constraints, it was also not possible to 
study the interaction of a proton-free ion beam with CD 
catchers. However, by comparing the data obtained from 
No.d^CD and CD^CD cases in Fig.O one can reason¬ 
ably assume a similar contributions from d(d, n)^He and 
d(p, n+p)^H reactions in the CD^CD case. 

Angularly resolved neutron spectra obtained by the 
six nToF detectors, as shown in Fig. [5l show the differ¬ 
ence in neutron generation between the cases. While the 
CD^C case produced low energy neutrons isotropically, 
similar to that obtained for pitcher only target shown in 
Fig. Efa), No-d^CD and CD^CD produced a very sig¬ 
nificant anisotropy in both neutron flux and maximum 
neutron energy, as shown in Fig.[5fb) and (c) respectively. 
Fig.E^c) shows the neutron spectra obtained along differ¬ 
ent angles for the case shown in Fig.[2fa), where the spec¬ 
trally integrated neutron flux over the neutron energy in 
the range 2.5-10 MeV obtained from the nToF detectors 
agrees well within the experimental errors with the off 
axis BDS and CR39 measurements shown in Fig. EJe). 
The highest neutron flux obtained in the experiment was 
for the CD^CD case, where, as shown in the Fig. [bfd), 
the nToF spectra agrees with that obtained from the BDS 
and activation diagnostics. The peak neutron flux along 
the beam axis, for neutron energy between 2.5-10 MeV, 
was close to 10^ n/sr. Indeed, the neutron flux can be 
significantly enhanced by optimising the neutron gener¬ 
ation with high yield catcher targets, such as lithium or 
beryllium, for which the reaction cross-section is an or¬ 
der of magnitude higher than for d(d,n)^He reaction and 
stays very high for ions with tens of MeV energy. 

The anisotropy ratios (which we define as the ratio 
of on-axis flux to off-axis flux) that have been observed 
in these experiments can be estimated for a beam-fusion 
scenario in terms of the projectile energy and the differen¬ 
tial cross-section of the nuclear reaction. For instance, in 
case of d(d,n)^He reaction, using the ion velocity in the 
centre-of-mass (c-o-m) frame, Vd^cm = (l/2)^£’c^/2mp, 
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FIG. 7. (a) Theoretical anisotropy ratio for d(d,n)^He reac¬ 
tions as a function of incident deuteron energy using Eq. [2 
(b) Energy of emitted neutrons in a beam-fusion scenario 
along 0 and 90 degrees to the incident ion beam axis as a 
function of incident deuteron energy, calculated using Eq. [T] 


one can obtain the emitted neutron velocity in the c- 
o-m frame from the energy-momentum conservation as 
Vn,cm = + Q)/4:mp, where is the incident 

deuteron energy in the laboratory (lab) frame, rUp is the 
mass of proton and Q is the Q-value of the reaction. 
Therefore, the neutron energy in the lab frame along a 
given neutron emission angle (6>, with respect to the in¬ 
cident ion beam axis) can be written as 


Eri — Ed /8 y/cos^ 0 2 6Q/Ed -1- cos 0 


1 2 


( 1 ) 


From this expression one can determine the neutron ve¬ 
locity in the lab frame Vn^ and thus the neutron emission 
angle in the c-o-m frame can be written as cos 0cm = 
{vn COS 0 — Vd)lvn^cm' The auisotropy is then determined 
from the differential cross-section in the lab frame, which 
is related to that in the c-o-m frame via, 

da (1 + -t- 2a cos 0cm)^^‘^ 

dCl 1 E a cos 0cm 



where, a = Vd,cml'^n,cm- An anisotropy ratio can thus be 
obtained at each deuteron energy by taking the ratio of 
the differential cross-section values that are calculated in 
this way at 0 and 90 degrees. This requires one to use the 
tabulated data for the differential cross-section in the c-o- 
m frame available in the Experimental Nuclear Reaction 
Database(EXFOR) ( 2 ^ . 

The results of this calculation are shown in Figl^^a), 
which indicates that the d(d,n)^He reactions could eas¬ 
ily produce the levels of anistropy shown in FigH] given 
that the deuteron energy spectrum shown in Figl6] ex¬ 
tends above 5 MeV. The anisotropy emission from this 
calculation assumes that the angular distribution of the 
incident ion beam is sufficiently narrower than that of 
the neutrons. Without the differential cross-sections for 
d(p, n+p)^H reactions in the c-o-m frame, it is hard to 
produce a similar graph for the d(p, n+p)^H reaction. 
However, even if one assumes an isotropic emission of 
neutrons for d(p, n+p)^H reaction in the c-o-m frame, 
the kinematic effects still lead to an anisotropy ratio of 
about 10 over a wide range of energies. 


In addition to the observed anisotropy in neutron flux, 
the energy distribution arising from the reactions, as can 
be seen from Fig. [5l is highly anisotropic due to nuclear 
reaction kinematics. According to Eq.[Tl one can approx¬ 
imate the energy of emitted neutrons along the beam for¬ 
ward direction (see Fig. EJb)) as En{0 = 0) ^ Ed while 
Ed > Q (3.3 MeV and -2.25 MeV for d(d,n)^He and 
d(p, n+p)^H reactions respectively), which is similar to 
the measured neutron energy by the nToF detectors as 
shown in Fig. [H 

In conclusion, we have demonstrated a strongly 
beamed (^ 70° FWHM), high flux (of the order of 10^ 
n/sr) source of fast neutrons based on beam-nuclear re¬ 
action employing high power laser driven protons and 
deuterons. The neutron flux in the beam, which was 
amongst the highest reported in the literature, was or¬ 
der of magnitude higher than that present elsewhere, as 
characterised spatially and spectrally by a suite of neu¬ 
tron diagnostics. Such a directed beam of fast neutrons is 
highly favorable not only for its direct applications, but 
also for its transport and moderation. With the possibil¬ 
ity of producing higher energy ion beams and higher flux, 
either by TNSA or via emerging ion acceleration mecha¬ 
nisms, and by employing higher yield neutrons converters 
(such as ^Li or ^Be), this approach can lead to develop¬ 
ment of an appealing neutron source for both established 
and innovative applications. 
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